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Abstract 

Cavia tschudii Fitzinger, 1867 is a wild guinea pig species living in South America that according to the 
analysis of mitochondrial genes is the closest wild form of the domestic guinea pig. To investigate the 
genetic divergence between the wild and domestic species of guinea pigs from a cytogenetic perspective, 
we characterized and compared the C, G and AgNOR banded karyotypes of molecularly identified Cavia 
tschudii and Cavia porcellus Linnaeus, 1758 specimens for the first time. Both species showed 64 chro- 
mosomes of similar morphology, although C. tschudii had four medium size submetacentric pairs that 
were not observed in the C. porcellus karyotype. Differences in the C bands size and the mean number of 
AgNOR bands between the karyotypes of the two species were detected. Most of the two species 
chromosomes showed total G band correspondence, suggesting that they ptobably represent latge 
syntenic blocks conserved over time. Pattial G band correspondence detected among the four submeta- 
centric chromosomes present only in the C. tschudii karyotype and their subtelocentric homologues in 
C. porcellus may be explained by the occurrence of four pericentric inversions that ptobably emerged and 
were fixed in the C. tschudii populations under domestication. The role of the chromosomal and genomic 
differences in the divergence of these two Cavia species is discussed. 
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Introduction 

Cavia tschudii Fitzinger, 1867 is a wild species of guinea pig (Rodentia, Caviidae) 
which inhabits northern Chile, southern Peru and Bolivia and northwestern Argentina 
(Weir 1974, Woods and Kilpatrick 2005). The domestic guinea pig Cavia porcellus 
Linnaeus, 1758 has a cosmopolitan distribution and is an experimental animal, pet, 
and even is consumed as food in countries of the Andean Altiplano (Tello 1972). 

There is a consensus that C. porcellus is a domestic form derived from one of 
the five currently recognized wild species of guinea pigs that inhabit South America 
(Woods and Kilpatrick 2005). The crosses between C. porcellus and C. fulgida Wa- 
gler, 1831 yielded offspring which behaved according to the Haldane's rule (Haldane 
1922), since females were fertile and males were sterile (Detletfsen 1914). By contrast, 
the crosses between C porcellus and C. aperea sensu Erxleben, 1777 (Pictet and Ferrero 
1951, Rood 1972) and between C porcellus and C. cutleri Tschudi, 1844 {sensu Ben- 
net, 1836) (Castle 1916) produced hybrids which were fertile in both sexes. C. aperea 
or C. tschudii have been repeatedly considered as the most probable ancestor of the 
domestic guinea pig. Later, molecular analyses of the mitochondrial cytochrome b and 
12S RNA genes clearly showed that the closest species to C. porcellus is C. tschudii and 
not the genetically related C. aperea (Spotorno et al. 2004, Dunnum and Salazar-Bravo 
2010). Based on these molecular results and on the analysis of mummified guinea pig 
remains found in archeological sites, Spotorno et al. (2007) suggested that the domes- 
tication of the wild guinea pig occurred in southern Peru-northern Chile. 

Considering that the karyotype provides useful characters in taxonomic and sys- 
tematic studies and that changes in the number and structure of chromosomes may 
contribute to speciation (King 1993, Searle 1993, Capanna and Redi 1994, Capanna 
and Castiglia 2004, Marques-Bonet and Navarro 2005, Faria and Navarro 2010), we 
describe and compare now for the first time the G, C and AgNOR banded karyotypes 
in molecularly identified specimens of the wild montane guinea pig C tschudii and the 
domestic guinea pig C. porcellus. Our objective is to discover the chromosomal and 
genomic differences between these two species of Cavia in relation to the divergence 
associated with the domestication process. 



Material and methods 

Skulls, skins and liver samples for DNA analysis of all the studied animals were pre- 
served in the collection of the Laboratorio de Citogenetica de Mami'feros, Facultad de 
Medicina, Universidad de Chile (LCM). We examined five Cavia tschudii specimens, 
four males (LCM 3199b, 3110, 3080, 3225) and one female (LCM 3232), collected 
in the locality of Molinos, Valle de Lluta, 18°23'S, 69°45'W Arica, I Region, Chile, 
and four Cavia porcellus animals, two males (LCM 2454, 3192) from the laboratory 
Pirbright breed, Instituto de Salud Piiblica, Santiago, Chile, and two females (LCM 
2479, 2489) from the Andean creole breed, Arica Agromarket, Arica, Chile. 
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Chromosomes were obtained from marrow cells using conventional in vivo col- 
chicine, hypotonic method, preceded by yeast injection to improve the mitotic index 
(Lee and Elder 1980). Metaphase cells were G-banded and C-banded by the meth- 
ods described by Chiarelli et al. (1972) and Sumner (1972), respectively. The active 
nucleolar organizing regions (NORs) were detected by the silver staining procedure 
(Sanchez-Rufas et al. 1982). At least 10 good-quality metaphases for each of the stain- 
ing methods per taxon were selected under a light microscope and digitally captured 
and stored. Chromosomes were counted, cut out and ordered by size and form using 
ADOBE PHOTOSHOP version 6.0. The centromeric indexes calculated by measur- 
ing the chromosomal arms in 12 metaphases of each species, allowed the classification 
of the chromosomes as metacentric, submetacentric, subtelocentric and telocentric 
(Levan et al. 1964). Chromosomes of both species were ordered in the groups de- 
fined by Fernandez and Spotorno (1968) for C. porcellus (groups A, B and C), add- 
ing a fourth group (group D) of submetacentric chromosomes for C. tschudii (Fig. 
1). Male and female G-banded karyotypes from each species were compared and the 
chromosomes were classified as having totally corresponding, partially corresponding 
or unique G band patterns (Spotorno 1977, Walker et al. 1979). The size and distribu- 
tion of the C and AgNOR bands were evaluated in six metaphases of C. tschudii and 
nine of C. porcellus. To determine the total number of active AgNOR in each of the 
species, AgNOR + sites were identified and counted in the chromosomes of 24 meta- 
phases per species. The statistical significance of the differences was estimated using a 
Chi squared test. 



Results 

Chromosome number, size and morphology 

For Cavia tschudii and Cavia porcellus we consistently found a 2n=64, FNa=100-102; 
the variation in the FNa of both species was due to the polymorphism of chromosome 
1 (Fig. 1). Cavia tschudii showed five pairs of submetacentric chromosomes (group D, 
Fig. la) of which four pairs (numbers 27, 28, 29 and 30) were not present in the C. por- 
cellus karyotype (Fig. lb). The X chromosome of C. tschudii was a large submetacentric 
similar to that of C. porcellus and the Y chromosome was a subtelocentric larger than 
that of C. porcellus (Fig. 1). 



G bands 

The comparison of C. tschudii and C. porcellus G-banded karyotypes revealed total 
correspondence for 25 of the 31 autosomal pairs and for the X chromosomes of both 
species (Figs. 2, 3, Table 1). The four submetacentric chromosomes present only in 
the C. tschudii karyotype showed partial G band correspondence with four C. porcel- 
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Figure I . Conventional stained karyotypes: a Cavia tschudii male b Cavia porcellus female. Chromo- 
somal pair N° 1 shows subtelocentric morphology in other individuals of both species (a and b insets). 
Bar = 5 pm. 
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Figure 2. G-banded karyotypes : a Cavia tschudii male b Cavia porcellus female. Chromosomes num- 
bered according to original karyotype descriptions (see Fig. 1). Bar = 5 pm. 



lus subtelocentric chromosomes (Fig. 6, Table 1). Only the Y chromosomes and two 
autosomal pairs (C. tschudii chromosomes 12, 13 and C. porcellus chromosomes 14, 
18) were unique of each species karyotype (Table 1). 



C bands 

The chromosomal distribution of the C bands was similar in the karyotypes of the 
two species, being located preferentially in the centromeres and the short arms of the 
chromosomes (Fig. 4). However, the amount of constitutive heterochromatin was ap- 
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Figure 3. Chromosomes of Cavia tschudii (Cts) and Cavia porcellus (Cpo) with total G band correspond- 
ence. Cts chromosomes are at the left and Cpo at the right of each chromosomal group. Note that the long 
arms of sub telocentric (St) and telocentric (T) forms of pair 1 show total G band correspondence (inset). 
Chromosomes numbered according to original karyotype descriptions (see Fig. 1). Bar = 5 um. 

Table I . Correspondence of Cavia tschudii and Cavia porcellus chromosomes according to their G band 
patterns 1 . 
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'Chromosome numbers are the one of each species karyotype (see Fig. 1); in the same column chromo- 
somes with total or partial G band correspondence. Cts = Cavia tschudii, Cpo = Cavia porcellus. 

preciably greater in C. tschudii than in C. porcellus, spreading over most of the short 
arms in several sub telocentric chromosomes (Fig. 4a). The X chromosomes of both 
karyotypes, equal in size, morphology and G bands (Figs 1-3), showed a C + band 
in the paracentromeric region of the short arm (Fig. 4). Both Y chromosomes were 
completely heterochromatic, being larger the Y chromosome of C. tschudii than the C. 
porcellus one (Fig. 4). 
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Figure 4. C-banded karyotypes: a Cavia tschudii male b Cavia porcellus male, showing heteromorphism 
for chromosome 1 . Most of the chromosomes of both species were tentatively identified according size 
and morphology. Bar = 5 um. 
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Figure 5. AgNOR-banded karyotypes: a Cavia tschudii male with four nucleolar chromosomal pairs 
(4, 11, 14 and 29) b Cavia porcellus male with five nucleolar chromosomal pairs (1, 3, 5, 7 and 1 1). The 
nucleolar chromosomes of both species were tentatively identified according to their size and morphology. 
Bar = 5 um. 

AgNOR bands 

Multiple AgNOR bands were detected in the karyotypes of both species, consistently 
located in the telomeres of several chromosomal pairs (Fig. 5). The analysis of some 
AgNOR banded metaphases per species indicated that the number of AgNOR bands 
was different between the two species and also among the individuals. Thus, the results 
showed that in C. porcellus the mean and maximum numbers of chromosomes with 
active NORs (5.76 and 9.0, respectively) were higher than those of C. tschudii (4.13 
and 7.0, respectively). Moreover, when we examined all the 3.072 chromosomes from 
48 metaphases of both species, each of them having 64 chromosomes, we found a total 
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Figure 6. Rough simulation of the changes associated with the occurrence of pericentric inversions in C. 
tschudii chromosomes. Submetacentric C. tschudii chromosomes (Cts, first column at the left) that origi- 
nate the subtelocentric C. porcellus chromosomes (Cpo, last column at the right): 1°) chromosomal break, 
2°) rotation of the cleaved segment, 3°) rejoining and sealing with the original segment. Chromosomes 
numbered according to original karyotype descriptions (see Fig. 1). Bar = 5 um. 



of 237 AgNOR + sites, 138 of them located in C. porcellus chromosomes and 99 in C. 
tschudii chromosomes. Accordingly, the number of chromosomes bearing active NOR 
was significantly higher in the C. porcellus karyotype than in the C. tschudii one (% 2 = 
6.956; p < 0.05; df= 1). 



Discussion 

Cavia tschudii and Cavia porcellus diploid numbers (2n=64), previously described with 
basic cytogenetic techniques (Ohno et al. 1961, Fernandez and Spotorno 1968, Dun- 
num and Salazar-Bravo 2006), were confirmed; nevertheless the fundamental number 
of autosomal arms (FNa = 100-102) were different to those reported before. For C. 
porcellus, Fernandez and Spotorno (1968) described an FNa = 96, while for C. tschudii, 
Dunnum and Salazar-Bravo (2006) found an individual in the Bolivian Altiplano with 
an FNa which ranged from 104 to 108. The FNa variability of guinea pig species may 
be due to polymorphisms for the presence of short arms in the chromosomes described 
as subtelocentric or telocentric in these species. The polymorphism for chromosome 1 
short arms detected previously for C. cobaya Pallas, 1766, a synonym of C. porcellus, 
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(Ohno et al. 1961, Schmid 1965, Zenzes et al. 1977) and for both species in this study, 
as well as the report of entirely heterochromatic short arms in the subtelocentric auto- 
somes of C.porcellus (Bianchi and Ayres 1971), give support to such hypothesis. 

The number and morphology of C. tschudii and C. porcellus chromosomes were 
similar to those reported for other subspecies and species of the genus Cavia. So, with 
the exception of C intermedia Cherem, Olimpio, Ximenez, 1999, and a population of 
C. magna Ximenez, 1980, having 2n=62 (Gava et al. 1998, Cherem et al. 1999, Gava 
et al. 2012), for all of the other taxa of the genus the same 2n=64 diploid number has 
been described, although with different numbers of autosomal arms. An FNa=124 was 
recorded for C. aperea pamparum (George et al. 1972); FNa=ll6 for C aperea aperea 
specimens from Pernambuco, Brasil (Maia 1984) and 114 for some individuals from 
the Bolivian lowlands (Dunnum and Salazar-Bravo 2006); FNa=124 for C. magna and 
C.fulgida (Pantaleao 1978) and 114 for C. nanaThomzs, 1917 (Cavia tschudii sodalis, 
1926) (Dunnum and Salazar-Bravo 2006). 

The analysis of the C bands showed that although they had a similar distribution in 
the chromosomes of the two species, they were smaller in size in the autosomes and in the 
Y chromosome of C. porcellus than in the C. tschudii ones, suggesting that a loss of het- 
erochromatin occurred during the domestication process. In accordance with this result, 
measurements of the genome sizes of 3 1 hystricognath rodent species (Gallardo et al. 2003) 
indicated that the genome of C tschudii (9. 1 pg) is larger than that of C. porcellus (8.2 ± 0.4 
pg), having the first species the largest genome size among the 30 diploid species analyzed. 

Five chromosomal pairs bearing NOR at the short arm telomeres were found by 
Zenzes et al. (1977) in the karyotype of the domestic Cavia, so being in agreement 
with our results. Using a double-staining procedure they could identify those chro- 
mosomes as numbers 1, 3, 9, 12 and 14 of the quinacrine banded stained karyotype. 
An accurate identification of the C. porcellus and C. tschudii nucleolar chromosomes 
described here would require the use of a similar double-staining procedure to allow 
the comparisons with other descriptions. 

The differences in the number of AgNOR bands found between the two Cavia 
species analyzed here and among the individuals in each of them, confirmed the ten- 
dency to variability in NOR expression usually described for mammals. It has been 
proposed that this variability would depend mainly on the specific metabolic demands 
of cells and individuals (Mikelsaar et al. 1977, Mayr et al. 1987, Sanchez et al. 1989, 
Suzuki et al. 1990, Bern'os et al. 1992, Zurita et al. 1997, Walker et al. 1999, Walker 
and Flores 2007). The greater number of AgNOR bands found in the C. porcellus 
karyotype than in the C. tschudii one would reveal a greater transcriptional activity 
of the ribosomal genes in the genome of the domestic form. One possible functional 
explanation of this result is that since C. porcellus has been selected for productive pur- 
poses, it would require higher rates of protein synthesis than the wild form. It should 
be mentioned that in a recent comparison of brain gene expression levels between four 
pairs of domestic/wild mammals, the largest differences were found between the do- 
mestic and wild guinea pigs, although C. aperea, genetically related to C. tschudii, was 
used as the wild guinea pig species (Albert et al. 2012). 
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Comparison of the G-banded karyotypes of the two Cavia species included in this 
study revealed that most of the autosomal pairs and the X chromosomes showed total 
G band correspondence, suggesting that these chromosomes constitute large syntenic 
blocks present in the common ancestor of both species and conserved over time. The 
differences in morphology and the partial G band correspondences detected between 
four chromosomal pairs of these two species, suggest that the four submetacentric chro- 
mosomes present only in the C. tschudii karyotype would have suffered pericentric inver- 
sions originating the four subtelocentric chromosomes of C. porcellus (Fig. 6, Table 1). 

Cavia porcellus would be the domestic successor of C. tschudii from which it would 
have originated more than 4000 and possibly 7000 years ago (Wing 1986) by a process 
of domestication and artificial selection in the C tschudii populations which inhabit 
southern Peru and northern Chile (Spotorno et al. 2007, Dunnum and Salazar-Bravo 
2010). While being domesticated, those populations must have been small in size 
and with only a few individuals participating as parents in the reproductive process, 
which over many generations would have produced high levels of endogamy. These 
characteristics would have facilitated the fixation of the pericentric inversions which 
must have emerged spontaneously and frequently in the populations. Specifically, the 
heterozygotes for the pericentric inversions would have decreased their fertility since 
their gametes would be unbalanced as a result of crossing-over in the inverted segment 
(Coyne et al. 1993, King 1993). As a consequence, gene flow between the original 
homozygotes and the homozygotes for the inversion would have been reduced, origi- 
nating genetic divergence between the two chromosomal forms. 

Nevertheless, it has been demonstrated recently that the fertility of the inversion car- 
riers is not always reduced (Muss and Schwanitz 2007). In some cases and depending on 
the size, genetic content, and chromosomal location of the inversion, the chromosomal 
inverted region pairs non-homologously with its normal partner forming a straight biva- 
lent which does not present any loop, so causing crossing-over suppression (Torgasheva 
and Borodin 2010). If that is the case, the absence or reduced recombination between 
the inverted and non inverted genomic regions in the Cavia pericentric inversions, would 
be the cause of genetic divergence accumulation and reduction of gene flow between the 
two chromosomal forms, as it was proposed as a general model of speciation by several 
authors (Noor et al. 2001, Navarro and Barton 2003, Hoffman and Rieseberg 2008). 

A critical assessment of reproductive isolation in crosses between Cavia species as 
previously reported, confronts the appropriate identification of specimens, the reliabil- 
ity of the taxonomy at the time, and the nature of the differences eventually found. For 
instance, crosses between C. porcellus and individuals from Arequipa, Peru identified 
as C cutleri Bennett, 1836, which correspond to the original description of C. cutleri 
based on a single specimen from lea, Peru (see Weir 1974), produced fertile offspring 
according to Castle (1916). By contrast, other wild specimens from lea, Peru that also 
received the name C cutleri by Tschudi in 1 849, were finally renamed as C. tschudii by 
Fitzinger in 1867 (see Weir 1974); the latter is now the usually accepted name for the 
wild montane guinea pig (Woods and Kilpatrick 2005). Therefore, C. cutleri Bennet, 
1838 is now considered a synonym of C porcellus (Woods and Kilpatrick 2005). If 
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Castle in fact crossed C. porcellus laboratory animals with C. cutleri Bennet (= C. porcel- 
lus) specimens, in reality he might be doing intraspecific crosses, and the fertility of the 
descendants would be an expected result. In any case, the assignation of the individuals 
from Arequipa to Cavia cutleri Bennet was not well documented in that study, since it 
was based only on the smaller body size of those individuals with respect to domestic 
C. porcellus (Castle 1916). Moreover, he did not indicated the mating times taken by 
the crosses, neither the number of pairs in which crosses were attempted, reporting 
only that a large number of descendants were obtained (n=107), as many as those ob- 
tained in crosses within each form (n=108). In sum, if the chromosomal and nucleolar 
differences we are reporting here in molecularly identified specimens of C. tschudii 
and C. porcellus were also found in other populations, we predict that their eventual 
hybrids will show some degree of genomic incompatibility. 

Reproduction of wild mammal species in captivity is a difficult and not always 
successful task. It is even more difficult to obtain descendants from crosses between 
different chromosomal races or species in the laboratory (Walker et al. 1984, 1999, 
Hauffe and Searle 1998, Castiglia and Capanna 2000, Franchini et al. 2008, Nunes 
et al. 2011). In crosses between phyllotine rodent species, we reported previously a 
decrease in the proportion of pairs with births and in the litter's size together with an 
increase in the time between mating and birth, compared to those registered for the in- 
traspecific crosses (Walker et al. 1984, 1999). Although we repeatedly tried to cross our 
specimens of C. tschudii with C. porcellus in our laboratory, we have had no success yet. 

If the analysis of the crosses realized between chromosomal races of Mus and Sorex 
rodents (Hauffe and Searle 1998, Castiglia and Capanna 2000, Franchini et al. 2008, 
Nunes et al. 201 1) would have only considered the number of descendants obtained, 
it would not have been possible to reach conclusions about the fertility level of those 
hybrids. To estimate fertility, specific reproductive aspects must be studied, such as 
the success obtained in crosses between the parental forms and some hybrid charac- 
ters, i.e.: the normality of their meiotic process, the histology of their gonads, the cell 
composition of their germinal line and the chromosomal constitution of the gametes 
that eventually they produce. In consequence, to evaluate the fertility level of eventual 
C. tschudii x C. porcellus hybrids, the reproductive characters just mentioned above 
must be analyzed in the descendants of crosses between individuals of the parental 
species taxonomically well identified. Specifically, the fertility of the heterozygotes for 
the pericentric inversions described here should be further investigated to evaluate the 
contribution of those chromosomal changes to the divergence of the two Cavia species. 
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